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MTI  Output  Detection  Statistics  for  Target 

Plus  Log-Normal  Clutter  \ 

i 

1.  INTRODUCTION  ! 

Generally,  analyses  of  the  detection  statistics  for  targets  immersed  in  clutter1  11  J 

do  not  include  the  effect  of  an  MTI  on  the  probabilities  of  detection  and  false  alarm. 

This  is  certainly  acceptable  for  targets  immersed  in  Rayleigh  clutter,  because  the 
clutter  output  from  the  MTI  is  still  Rayleigh,  except  with  a different  variance.  How- 
ever, it  is  unacceptable  for  log-normal  clutter.  That  is,  if  the  input  to  the  MTI  is 
log-normally  distributed  the  residue  at  the  MTI  output  is  not  log-normally  distributed. 

In  this  report  we  will  therefore  derive  the  detection  statistics  for  the  output  of  a two- 
pulse  MTI  canceller,  when  the  input  to  the  MTI  consists  of  a target  plus  log-normally 
distributed  clutter.  Because  of  the  mathematical  complexity  we  shall  limit  our 
analyses  to  the  two-pulse  canceller,  but  the  generalization  to  the  N -pulse  canceller 
is  conceptually  straightforward. 

(Received  for  publication  17  August  1978) 
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As  a model  for  our  canceller  we  shall  consider  the  quadrature  channel  system 
shown  in  Figure  1.  For  purposes  of  review  we  will  first  derive  the  residue  proba- 
bility density  function  for  the  case  when  the  clutter  amplitude  is  Rayleigh.  We  will 
then  consider  the  more  difficult  case  of  log-normally  distributed  clutter,  and  will 
derive  single-hit  detection  statistics  (that  is,  probabilities  of  false  alarm  and  detec- 
tion) for  the  cases  of  a constant  target  and  a Rayleigh  distributed  target. 


IN  PHASE 
COMPONENT 
(TARGET* 
CLUTTER) 


QUADRATURE 
COMPONENT 
(TARGET  + 
CLUTTER) 


Figure  1.  Quadrature  Channel  Canceller 


2.  ANALYTICAL  PRELIMINARIES 

The  return  at  time  t from  a target  immersed  in  clutter  can  be  written  for 
narrow-band  signals,  as 

z(t)  = s exp  ^j(wdt  + 0)^j  + C exp  |^j(coo t + ^)j  (1) 

where  s and  <4  are  the  amplitude  and  phase  of  the  voltage  produced  by  the  target, 
is  the  doppler  frequency  of  the  target,  C and  \j/  are  the  amplitude  and  phase  of  the 
voltage  produced  by  the  clutter  and  u is  the  doppler  frequency  of  the  clutter.  At 
a later  time  t + T the  received  signal  can  be  written  as 


v 


6 


(2) 


z(t+  T)  = s'  exp  j^ju>d(t+  T)  + j0  'J  + C'  exp  ^ja)Q(t+  T)  + jV'  ' j 

•ff 

where  s'  and  <f>'  are  the  target  voltage  amplitude  and  phase  at  time  t+T.  If  we 
assume  that  the  target's  amplitude  and  phase  do  not  vary  from  pulse-to-pulse,  so 
that  s'  = s and  <p'  = tp,  we  then  can  readily  see  that  the  output  of  a two-pulse  can- 
celler is 


w T 

A = z(t)  - z(t+  T)  = 2 s sin  — 5 exp 


+ (C  e^  - C'  e 


j*'» 


1(0  + wdt  +-£-  T - j'j 


exp  (j  «Qt)  . 


(3) 


We  now  assume  that  T = set  </>  + co^t  = 2njr,  and  also  assume  that  there  is  no 

clutter  motion,  so  that  u)Q  = 0.  We  can  then  write  the  in-phase  and  quadrature  com- 
ponents of  A = r exp  (j  9 ) as 

r cos  9 = 2 s + C cos  0/  - C'  cos  ' 

r sin  9 = C sin  \p  - C'  sin  ^ ' . 

Our  task  now  is  to  determine  the  probability  densities  of  the  residue  amplitude,  r, 
and  phase  9 , in  terms  of  the  probability  distributions  of  the  target  and  clutter.  We 
will  first  consider  the  rather  simple  case  when  the  in-phase  and  quadrature  compo- 
nents of  the  clutter  are  gaussian  random  variables,  and  then  proceed  to  the  more 
complex  case  of  log-normall'  distributed  clutter. 


(4) 

(5) 


3.  RAYLEIGH  CLUTTER 

Let  us  study  the  case  when  the  clutter  voltage  amplitude  is  a Rayleigh  random 
variable.  This  implies  that  the  in-phase  and  quadrature  components  of  the  clutter 
are  joint-gaussian  random  variables.  Let  us  rewrite  Eqs.  (4)  and  (5)  as 


x = r cos  9 = 2 s 

+ c - c'  . 

(6) 

X X 

y = r sin  9 = C 

- c'  . 

(7) 

J y 

y 

Then  if  C , C , C',  C'  are  zero-mean  gaussian  random  variables  we  can  write  their 

* y * y 

joint  probability  density  as 


Also  C'  s C(t+T),  i//'  =iMt+T). 
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t 

> 

I 


p (C  , c , c 
*c  x y 


* cf  *c;2  - !p,CiC'  - 2ppcyc;]} 


C 2 . C..2 


(8) 


where 

P 


. <c.o  - <v;> 
1 <7»7>  <v> 


»,2  ■ <^>  * <s2>  - <c">  ■ <c;>  • 

In  writing  Eq.  (8)  we  have  assumed  that  the  in-phase  and  quadrature  components  are 
uncorrelated.  The  justification  for  this  assumption  is  given  in  Beckmann.  The 
probability  density  of  (x,  y)  keeping  s fixed  is  then  given  by 

oo  °o 

P(x.y|s)  = dC;  Jx  dC'  pc  [cx=x-2s+c;.  Cy  = y + C^.  C'.  C'  ] 

Upon  substituting  Eq.  (8)  into  Eq.  (9),  and  performing  the  integrations,  we  get 


. (9) 


p(x,y  I s)  = ^ 


i (x  - 2 s)2  + .y2 
- exp 

4 7ja^  (1-p)  ( 4(1  - Pj'dj 


(10) 


We  now  recall  that  x = r cos  0 , y = r sin  9 . Therefore  p(r,  0 I s)  is 
p(r,  0 | s)  =- 


1 

_( 

^r2  + 4 s2  - 4 r s cos  0^ 

exp 

4 TTOj  (1  -Pj) 

V 4a1*<1-Pl)  J 

1 

(11) 


If  we  integrate  Eq.  (11)  on  0 from  -n  to  v we  then  get 
p(r  | s) 


exp 


2o  x <1  -pj) 


2 p 2 

r + 4 s 
40^(1  - p1> 


'•  (7^) 


(12) 


T.  Beckmann,  P.  (I96  0 Probability  in  Communication  Engineering,  Harcourt, 
Brace,  and  Jovanovich,  New  York. 

8.  Davenport,  W.  , and  Root,  W.  (1958)  Random  Signals  and  Noise,  McGraw-Hill 
New  York. 

Note  that  s is  the  signal  level  assuming  the  target  return  is  solely  at  frequency  u>  . 
If,  as  is  often  done,  the  target  signal  is  uniformly  averaged  over  radial  velocities, 

1/2  22. 

then  we  must  replace  s by  s / ( 2)  irt,Eq.  (12),  and  by  <JQ  / 2 in  Eq.  (14).  The 

factor  of  2 comes  from  averaging  sin  (u)dT/2). 


8 


Equation  (12)  is  the  probability  density  of  the  residue  amplitude  for  a constant 
target  amplitude.  If  the  target  amplitude  fluctuates  and  these  fluctuations  obey  a 
a probability  density 


we  find  that 


OO 

p(r)  = J p(r  | s)  p(s)  ds 
o 


r exp  , T 

4[CT1  (1-p^*  20o^] 

2[a1a(i -Pl)  +2  oj*] 


(14) 


- 1 9 

The  quantity  (1  -p^)  has  been  defined  by  Nathanson  as  the  clutter  improvement 
factor  I.  Consequently,  for  Rayleigh  clutter,  we  see  from  Eqs.  (12)  and  (14)  that 
the  effect  of  the  two-pulse  canceller  is  to  reduce  the  clutter  variance  by  the  clutter 
improvement  factor  I =[  1 - p j(T)]  The  results  in  Eqs.  ( 12)  and  (14)  are  not 
unexpected,  and  have  been  presented  only  for  completeness. 
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4.  LOG-NORMALLY  DISTRIBUTED  CLUTTER 

In  order  to  consider  the  case  when  the  clutter  is  log-normally  distributed  let 
us  write  C(t)  = exp(x)  and  c'  = C(t+T)  = exp(x*),  so  that  Eqs.  (4)  and  (5)  become 

/ 

r cos  9 = 2s  + e^  cos  - e^  cos  \l/‘  , (15) 

/ 

r sin  9 = e^  sin  i p - sin  i//  , (16) 

We  now  assume  that  , t//  and  \j/'  are  jointly  distributed  gaussian  random 

variables.  That  is,  if  we  define  yj  = X»  yg  = VG  yg  = X and  y^  = \p'  their  joint 
probability  density  is  given  by^ 


9.  Nathanson,  F.  (1969)  Radar  Design  Principles,  McGraw-Hill,  New  York. 

$ 

Note,  that  this  differs  considerably  from  Rayleigh  clutter.  In  that  case 
C = exp(x)  cos  i//,  C = exp(x)  and  i//  and  so  on,  were  gaussian  random  variables. 

x y 

Here,  it  is  xand  ^ themselves  which  are  gaussian. 
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p(>v  y2*  y 3*  y4} 


4 4 


exp 


— A-rS  U 

2|  A|  n=  1 m=  1 


i ^n  " <yn>)  (yr 


-TTT 


( 2 77  ) 


(17) 


where 


|A|nm  is  the  cofactor  of  the  element  Anm  in  the  determinant  | A | of  the  covariance 
matrix,  and  ( > denotes  an  ensemble  average.  The  form  of  the  probability  distribu- 
tion of  the  residue  amplitude,  r,  and  phase  0 , for  the  general  case  when  \,  x*  ^ 
and  \Js ' satisfy  Eq.  (17)  is  exceedingly  complex.  The  form  of  these  results  can  be 
simplified  considerably  if  we  make  the  following  assumptions: 

(1)  The  mean  phase  angle  of  the  clutter  return  is  zero.  That  is, 

<*>=  <^'>  = 0 , 

(2)  The  phase  \p  and  log-amplitude  x are  uncorrelated.  That  is, 

«X'  <X>>  = <^-  <(!/'>)>  <<^>  = 0, 

(3)  The  correlation  function  p(T)  is  the  same  for  both  the 
log-amplitude  and  phase  fluctuations.  That  is, 

<<X-<X>)  (x'  - <x'>  )>  =cr  2p(T), 

A 

<< P = O^2  p(T) . 


In  those  cases  where  the  aforementioned  assumptions  may  be  made,  Eq.  (17)  re- 
duces to 


p(x.  x'»V>  'P>)  = P,(X,  x'>  P2<^(  'p') 


(18) 


where 


exp 


Pj(X.  x')  = 


_ [(X-  <X>>2  - 2p(x-<x~>>  (x'-<X^  > + (x'-<\*>  )2] 


2ax  (1-p  ) 


2 770^  <1  -pV/2 


(19) 
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exp 


L [ ^2-2p  w' + I//'2  j 


p2^.^')  1 2V(i-^2> 


217(7^(1  -p2)1'2 


(20) 


and  the  ranges  of  the  phase  angles  are  -°°  s \p  s oo  and  -"osij/'s*  Ordinarily, 
we  limit  the  phase  to  the  intervals  -77  ^ 4/  ' s 77,  -77  ^ s 77,  rather  than  to  the 
range  to  °°.  If  we  do  this,  we  have  shown  in  Appendix  A that  for  p r 1,  Eq.  (20) 
becomes 


exp 


pJ<P,  41  ) ~ 


4(1  -p)ff. 


4 773/2cr^(l  -p) 


(21) 


where  now  -77  £ if/  s 77,  -77  £ !^/  =»  77. 

Now  that  p(x»  X*.  ty')  is  known,  the  next  step  is  to  determine  p(r,  9 , x*.  41  ' I s). 
In  order  to  do  this  we  first  define  x*  = In  v and  set  <x>  = <X<>  = (?n  7;  . Also 
we  use  Eqs.  (15)  and  (16)  to  express  x and  4^  in  terms  of  r,  9 , v and  i p'  for  a given 
signal  level  s.  We  get 


X = -|  ^n  g , 


g=[r  +4s  + v + 2rv  cos  (9  - 1//  *)  - 4 s r cos  9 - 4 s v cos  ^ '] 


4/  = h = tan 


1 ( r sin  9 + v sin  l// ) 

( r cos  0 - S s + v cos  ip ' f 


(22) 

(23) 

(24) 


and 


dX  d*  dX'  d*'  = rdfdMv  W 


(25) 


If  we  use  the  results  of  Eqs.  (22)  to  (25)  in  Eqs.  (18)  to  (20)  we  find,  assuming  that 
the  correlation  coefficient  p1  T)  r 1, 

p(r,  0,  v,  t//|s)  = r exp  - 5-  | [ln&- - Zn— — Y 

( 4(1  - p)  a 2 L \ Vo  Vo  ) 


+ 2(1  - p)  i!n 


(26) 
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Finally,  we  obtain  the  probability  density  for  the  residue,  r (holding  s fixed) 


7T  7T  oc 


p(r|s)=  g »w-  r w J di//  j d6  J 

2<2t7)5/2<1 -p)a2a.  ^ ^ 


/ , r / 1/2  \ 2 1/2' 

< jr  Mn&- J «■  2<1  -P>  «n-S- — 

{ 4(1- p)  a (_\  ^ o ^°/  **  o **  o _ 


4>  U-P) 


(h-V/  ) 


In  order  to  perform  the  integrations  in  Eq.  (27)  we  first  realize  that  because 
2 - 1 2 - 1 

1 - p - 0 the  quantities  [4a^  (1  - p)]  and  [4  V (1  -p)]  in  the  exponential  are  very 

large.  Consequently  the  integrand  is  extremely  small  unless  ( in  ^ — - in~—  ) 

,2  V 7<  o I o / 

and  (h  - \p  ) are  quite  small.  This  condition  occur^  only  for  v > > r and  v > > s,  so 

that  we  may  expand  (h  -i//)2  and  (in  - in  — — ) in  a Taylor  series  in  the  small 

V o ho 

quantities  r/v  and  s/v.  If  we  ignore  higher  order  terms  in  these  small  quantities 
we  find 

h -xj/'  - — sin  (9  -ty')  + sin  i p'  , (28) 


in  — - in ~ — cos  (6  - \f/  ) cos  , 

77  7)  v v 

' o o 


g r v . (30) 

If  we  now  use  Eqs.  (28)  to  (30)  in  Eq.  (27),  and  further  assume,  for  mathematical 
convenience,  that 

a 2 = a , 2 = a2  (31) 

X ^ 


we  obtain  (again  ignoring  higher  order  terms  in  r/v,  s/v) 


p(r  | s)  e/g  r / d6  J ' dtf/ ' J exp 

2(277)5/‘i  (1  - p)CT3  -77  -77  O V3 


( 4 cr  ( 1 - p)  L 


cos  0 + 


211  ,n!t:]| 


4 


t 


The  i//'  integration  in  Eq.  (32)  is  trivial,  and  the  integral  on  8 is  immediately 

recognized  as  the  definition  of  the  modified  Bessel  function  IQ(.  . . ).  Consequently, 

if  we  do  these  integrations,  and  then  make  the  transformation  of  variables 
2 2 

y = T)  /v  we  obtain 


p(  r | s)  = 


4(2n)1 1 2 a3 1]  2 ( 1 - p) 
o 


exp 


(33) 


Finally,  we  define  new  normalized  residue  and  signal  variables 
R s L. 


2T)  0(1  -p) 

o 


TTT  ’ 


S = 


2r]  o<l-p) 

o 


TTT  * 


(34) 

(35) 


In  terms  of  these  new  variables  Eq.  (33)  assumes  the  simpler  form 

p(R I S)  = — J dy  I (4  S Ry)  exp  j-  (R2  + 4 S2)  y - ^ n J > . 

(27T)i/^o  n U ( 8 O*2  ) 


(36) 


0 s R S oo 


Equation  (36)  is  the  normalized  two-pulse  canceller  residue-amplitude  distribution 
for  a constant  target  immersed  in  log-normal  clutter. 

We  can  also  obtain  the  residue  probability  density  for  a fluctuating  target.  For 

% 

example,  suppose  the  target  fluctuations  have  a probability  density  (0  s s < °°) 


Jj  exp  (-^fr) 

r\  ' r\  ’ 


(37) 


We  emphasize  that  S is  the  signal  level  assuming  the  target  return  is  at  uq  only. 

u 2 

If  the  target  were  uniformly  averaged  over  radial  velocities  we  would  replace  S 
by  S*  12  in  Eq.  (36).  Likewise  in  Eqs.  (39)  and  (40),  we  would  replace  aQ2  by 
Cq2I2. 
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Then  in  terms  of  S 


4(1 -pit)  2 a2 

p^(S)  = 2 — S exp 


2(1  -p)  t)  2 a2  S2 


(38) 


Upon  multiplying  Eq.  (36)  by  Eq.  (38)  and  then  integrating  on  S from  0 to  00  we  obtain 


iiriW 


/ u \ H ( J F l 1 + 4 y y 

P(K)  = — J/2  J UL- 


(277)"  “a 


1 + 4 y y 


(39) 


where 


y = 


2 — • 

2(1  - p)  7]  a 
' o 


(40) 


Equation  (39)  is  the  probability  density  of  the  normalized  two-pulse  canceller- 

residue  amplitude  for  a fluctuating  target  immersed  in  log-normal  clutter. 

For  mathematical  convenience  we  next  assume  that  the  residue  voltage  is  passed 

through  a square-law  detector.  If  we  define  the  output  of  the  square-law  detector 
2 2 

as  Z = R and  also  define  W = S then,  for  a constant  target,  the  square-law  de- 
tector output  amplitude  satisfies  a probability  density  function 

oo 

p(Z  | W)  = X-TTr  J dy  I [4y(ZW)1/2]  exp  j-  (Z+4W)y  - (41) 

2a (277)  1 ' ^ 0 0 I 


For  a fluctuating  target,  the  square-law  detector  output  Z has  a probability  density 

2 

~ AT 

1 


p(Z,  y)  = 


2(277)  i/2a  "o 


J dyeXp  ll+TTV 


&n  y 

8cr2 


1 + 4 y y 


(42) 


Before  presenting  numerical  results,  it  is  informative  to  express  Z,  W and  y 

in  terms  of  the  clutter  improvement  factor,  I.  In  Appendix  B we  show  that  1 is 

2 

related  to  the  variance,  a , and  correlation  coefficient  p,  of  the  log-amplitude 
fluctuation  as 


I=2CA~  -7 


1 


2a  (1  -p) 


(43) 
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where  is  the  clutter  attenuation  factor,  as  defined  by  Nathanson  . Consequently, 

we  can  write  Z = R2  = ( r 2 / 2 n 2)  1,  W = <s2 / 2rj  2)  1 and  y = (a  In  )2  I,  where  n 

o o o o o 

is  the  median  clutter  amplitude  at  the  input  to  the  two-pulse  canceller.  Thus  y 
represents  the  improved  signal  to  median- clutter  ratio  for  a fluctuating  target  at 
the  center  doppler  frequency  = tf/T  . 

The  integral  in  Eq.  (42)  has  been  evaluated  numerically  for  a number  of  differ- 
ent values  of  y,  and  the  results  are  presented  in  Figure  2.  We  observe  that  the 
target  free  (y  = 0)  residue  probability  density  function  has  a very  different  behavior 
than  that  for  the  case  of  large  signal-to-clutter  ratios  (y>  > 1).  This  will  be  clearly 
evident  in  the  next  section,  where  we  evaluate  p(Z,  y)  in  Eq.  (42)  approximately. 
Similar  results  for  the  case  of  a constant  target  (that  is,  from  Eq.  (41))  are  shown 
in  Figure  3. 

It  is  also  important  to  obtain  results  for  the  probabilities  of  false  alarm  and 
detection.  The  probability  of  false  alarm  is  defined  as 
00  ) 

Pf  = J dZ  p(Z,  y = 0)  , (44) 


and  for  a fluctuating  target  is  plotted  in  Figure  4 for  the  case  when  CT=  0.  707. 


Figure  2.  p(Z,  y)  for 


The  probability  of  detection  is  defined  via 

oo 

Pd  = £ dZ  p(Z,  y)  . (45) 

o 

where  p(Z,  y)  is  given  by  Eq.  (42).  Results  for  are  plotted  in  Figure  5,  for  the 
case  when  a = 0.707. 


5.  APPROXIMATE  FORMULAE  FOR  A FLUCTUATING  TARGET 


In  this  section  we  will  present  approximate  formulae  for  p (Z,  y)  P^.  and  P^. 
When  the  target  is  absent  (y  = 0)  it  is  possible  to  evaluate  Eq.  (42)  for  large  values 
of  the  normalized  residue,  Z,  by  the  method  of  steepest  descent.  The  result  is 


p(Z. 


y = 0)  = 


e'1 

~2cT 


exp 


(46) 


2 

where  v = 4 a and  ( = i<+  fn  i/Z.  A comparison  of  results  obtained  from  Eq.  (46) 
with  the  exact  computer  evaluation  of  Eq.  (42)  is  given  in  Table  1,  for  the  case  when 
a = 0.  707.  Observe  that  Z > 1,  Eq.  (46)  is  a very  good  approximation  to  Eq.  (42). 
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Table  1.  Comparison  of  the  Approximate  Result  in  Eq.  (46)  With  the 
Exact  Evaluation  of  Eq.  (42) 


z 

p(z,  y = 0)  — from  Eq.  (42) 

p (z,  y = 0)  - from  Eq.  (46) 

0.  1 

1.  2512 

0.  9131 

0.  3 

0.  6456 

0.  6667 

1 

0.  2134 

0.  2197 

3 

0.  0517 

0.  0528 

10 

0.  006669 

0.  006760 

30 

0.  000634 

0.  000639 

100 

0.  0000275 

0. 0000276 

300 

9. 294  X 10'7 

9.  294  X 10'7 

For  Z -+<*>'  Eq.  (46)  can  be  approximated  by  the  simpler  form 

1.1 


p(Z,  y = 0) 


e 1 Un  vZ) 


- + 


(rZ) 


1 + 1/r 


exp 


(47) 


Thus,  even  for  Z-*«>  the  residue  pdf  is  clearly  not  log-normal,  although  it  has  a 
relatively  slow  dropoff  rate  which  is  a characteristic  of  the  log-normal. 

Next  we  consider  the  pdf  when  the  output  signal-to-clutter  ratio  is  large.  When 
y » 1 we  can  again  obtain  an  analytic  approximation  by  ignoring  unity  in  compari- 
son with  4 y y in  the  integral  in  Eq.  (42).  The  result  is 


p(Z,  y >>  1)  tr  ^ 1 


exp 


(48) 


Results  obtained  using  Eq.  (48)  are  compared  with  the  exact  evaluation  of  the 
integral  in  Eq.  (42)  in  Table  2.  Note  that  the  agreement  is  quite  good  provided 
Z«  4 y2. 

A comparison  of  Eqs.  (46)  and  (48)  clearly  shows  that  the  nature  of  the  residue 
probability  density  function  is  quite  different  for  large  values  of  target-signal  to 
clutter  ratio  than  for  the  case  when  the  target  is  absent;  in  the  former  case  the  pdf 
is  exponential  whereas  in  the  latter  case,  it  is  a variant  of  the  log-normal. 


We  have  also  obtained  simplified  expressions  for  the  probabilities  of  false 
alarm  and  detection.  For  large  thresholds  (Zq  >>  1)  we  find 


r 


■ 


where  /f  = 0n  (1/  Zq>  and  ^ = 4«^,  Comparisons  of  the  results  obtained  using  Eq.  (49) 
with  exact  computer  calculations  of  Pf  are  given  in  Table  3.  We  observe  that  the 
agreement  is  excellent. 


Table  2.  Comparison  of  the  Approximate  Result  in  Eq.  (48) 
With  tile  Exact  Evaluation  of  Eq.  (42) 


> 

= 10 

7 

= 100 

z 

p(z,  y) 
exact 

p(Z,  y ) from 
Eq.  (48) 

P(Z,  y) 
exact 

p(Z,  y ) from 
Eq.  (48) 

0.  1 

0.  0237 

0.  0250 

0.  00248 

0.  00250 

1 

0.0231 

0.  0243 

0.  00248 

0.  00249 

10 

0.  0186 

0.  0195 

0.  00242 

0.  00244 

100 

0.  00222 

0.  00205 

0.  00194 

0.  00195 

300 

0.  0000245 

0.  0000138 

0.  00118 

0.  00118 

1000 

0.  000207 

0.  000205 

Table  3.  Comparison  of  the  Exact  Values  of  Pf  With  the 
Approximate  Results  Obtained  Using  Eq.  49)  for 
(7=0.  707 


Zo 

Pf  (exact) 

Pf  (approximate  - from  Eq.  (49) 

314 

740 

1591 

1.  0001  X 10~4 

1.  0072  X 10~5 

1.  0013  X 10'6 

1.  0070  X 10‘4 

1.0120  X 10'5 

1.  0044  X lo"6 

When  Zq  is  sufficiently  large  and  v - 4(7  is  of  order  unity,  we  can  approximate 
Eq.  (49)  further  as 


_1  _ 1 

P{*(v  ZJ'111' Un  Zq)  p 1 exp 


in 


(50) 


It  is  interesting  to  speculate  what  results  we  would  have  obtained  if  we  had 
simply  assumed  that  the  residue  clutter  distribution  were  log-normal,  with  the  same 
variance  as  the  input  clutter  distribution,  but  with  the  median  level  squared,  tj  ^ , 
divided  by  the  clutter  attenuation  factor  C'A  = 1/2.  Then  upon  using  Eq.  (30)  in 
Reference  4,  witht/^2  replaced  by  2 T)  2/I,  we  would  find 
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p 


o 

f 


1/2 


exp 


IHT2TT 

o 


In 


(51) 


A comparison  of  the  result  in  Eq.  (51)  with  the  exact  results  for  Pj.  indicates  that 
this  is  a poor  approximation,  and  leads  to  threshold  levels  which  are  incorrect  by 
a factor  of  2 or  more. 

Finally,  for  large  values  of  the  equivalent  signal-to-clutter  ratio,  y , we  can 
derive  a suitable  approximation  for  P^.  This  is 


where  is  the  threshold  required  to  give  a specified  value  of  P^.  A comparison 
of  results  obtained  using  Eq.  (52)  with  exact  computer  calculations  of  is  given 
in  Table  4.  We  see  that  Eq.  (52)  is  quite  accurate  whenever  y is  sufficiently  large 
that  P^  > 0.  05.  Because  this  is  usually  the  regime  of  interest,  it  is  evident  that 
Eq.  (52)  is  a useful  approximation. 

Although  we  have  plotted  numerical  results  only  for  a = 0.  707,  it  is  clear  that 
Eq.  (49)  and  Eq.  (52)  can  be  used  to  obtain  results  for  other  values  of  a.  Some 
typical  results  are  shown  in  Figures  6 to  9.  Observe  that  the  detection  statistics 
depend  strongly  on  CT. 


Table  4.  Comparison  of  Exact  Computer  Calculations  of  P , With  the 
Approximation  in  Eq.  (52)  fora  = 0.707 


pf= 

io'4 

Pf=  10'5 

Pf=  10'6 

Pd 

m 

■H 

jhh 

Pd 

(Eq.  52) 

(Eq.  52) 

10  0.00103 

30  0.0773 

0. 00039 

0.  0730 

0.  0000203 

0.  00248 

9. 24  X 10'9 

0.  00210 

6. 52  X 10"6 

1.  74  X 10'6 

100  0.458 

0.456 

0.  1589 

0.  1572 

0.  0192 

0.  0187 

300  0.7698 

0.  7698 

0.  540 

0.  540 

0.  266 

0.  266 

i 
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Figure  6.  Probability  of 
F alse  Alarm  as  Computed 
From  Eq.  (49) 


Figure  7.  Probability  of 
Detection  for  Pf  = 10“  5 as 
Computed  From  Eq.  (52) 


6.  DISCUSSIONS 


I 

r 

I 

1 

* 


It  is  rather  interesting  to  compare  the  detection  statistics  for  log-normally 

distributed  clutter  with  those  for  Rayleigh  clutter.  Suppose  we  assume  that  the 

median  clutter  level  is  the  same  for  both  the  Rayleigh  and  log-normal  clutter  returns, 

and  we  also  use  the  same  definition  of  y for  the  Rayleigh  distribution  as  for  the  log- 

2 2 

normal  case.  That  is,  we  again  define  y = (Oq  /?)  q ) I,  where  I is  the  clutter  im- 
provement factor.  If  we  assume  that  the  residue  voltage  in  Eq.  (14)  is  passed 
through  a square  law  detector,  define  R^  = r I/27)q  , and  then  integrate  the  result 
to  obtain  and  P^.,  we  find  that  for  a fluctuating  target  immersed  in  Rayleigh  dis- 
tributed clutter 


Pd  = exp 


1 + 


ta  Pf  I 

2?v'2  ) 


(53) 


For  Rayleigh-distributed  clutter  a j,  and  7)  are  related  via  7] 
Consequently  Eq.  (53)  becomes 


exp 


in  P 

1 + 5.  545 


2.  7725  CTj2. 


(54) 


The  result  in  Eq.  (54)  is  plotted  in  Figure  10  for  a probability  of  false  alarm, 
Pj.  of  10  6 . Also  plotted  on  the  same  figure  is  the  probability  of  detection  for  the 

case  when  the  clutter  is  log-normally  distributed,  with  the  same  median  value,  7)  , 

2 ° 
as  for  the  Rayleigh  clutter,  and  variance  a = 0.  5.  We  note  that,  because  the 

signal-free  residue  pdf  and  consequently  the  threshold,  is  very  different  there  is  a 

considerable  difference  in  the  probability  of  detection  in  log-normal  as  opposed  to 

Rayleigh  clutter. 

Our  discussion  here  has  been  limited  strictly  to  the  case  of  a single  hit  (that 
is,  a single  output  from  a'two-pulse  canceller).  If  the  target  and  clutter  outputs 
remain  nearly  perfectly  correlated  from  pulse-to-pulse  then  the  M-hit  probabilities 
of  false  alarm  and  detection  are  the  same  as  those  for  one  hit. 

We  have  also  limited  our  discussion  to  the  two-pulse  canceller.  The  extension 
of  our  analysis  to  an  N-pulse  canceller  is  straightforward  but  extremely  laborious, 
and  will  not  be  pursued  here. 
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7.  CONCLUSIONS 

We  have  derived  the  output  probability  densities  and  detection  statistics  for  a 

two-pulse  canceller  when  the  input  to  the  canceller  consists  of  a target  immersed 

in  log-normally  distributed  clutter.  We  have  shown  that,  even  with  the  target 

absent,  the  output  residue  is  not  log-normally  distributed  but  has  a probability  den- 

2 2 

sity  which  can  be  approximated  by  Eq.  (46)  for  large  values  of  Z = r I/2t)  q . This 
distribution,  however,  does  have  a relatively  slow  dropoff  with  the  consequence  that 
a relatively  high  normalized  thresholds  Zq,  is  required  in  order  to  achieve  a given 
probability  of  false  alarm.  This  threshold  is  much  higher  for  log-normal  clutter 
than  it  is  for  Rayleigh  clutter. 

Finally,  we  again  emphasize  that  our  definition  of  signal  power  level  always 

assumes  that  the  return  is  strictly  at  the  center  frequency  U).  = ff/T.  For  an  average 

2 2 a 

over  radial  velocities  we  would  replace  s by  s /2  in  Eq.  (12),  Eq.  (36),  and  so  on. 

2 2 

Similarly  in  Eqs.  (14),  (39),  (40),  and  so  on,  we  would  replace  <JQ  by  Oq  / 2 to 
include  the  effect  of  the  velocity  average. 
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Appendix  A 

Derivation  of  Equation  (21) 


We  consider  the  probability  density  giver,  hv  Eq.  (20)  which  is  defined  for  p 
and  i p'  ranging  from  to  °°,  and  realize  that  if  we  limit  the  range  of  the  phase 
from  -77  to  77  then  the  probability  of  the  phase  lying  in  dp  d p ' about  (p,  p'),  where 
1 p,  p'  are  now  limited  to  the  range  -77  to  77,  is  equal  to  the  probabilities  of  (p , p') 
in  Eq.  (20)  lying  in  d xp  dp'  about  (p  + 2 n77,  p'  + 2 m!7)  for  all  possible  integers 
n and  m.  That  is 


p(p,  p’)  = 


E E 


277 (1-p  ) n=-°°  m=-°° 


+ 2(l-p){p  + 2n7?)  (p‘  + 2 m 


exp 


1 

5 ~~T 

2o^(\-p*) 


*>]  • 


\(p  + 2 n 77  - p'  - 2 m77r 


(Al) 


where  p and  p’  now  lie  between  -77  and  77.  Now  if  the  correlation  from  pulse  to  pulse 

is  high,  as  is  necessary  for  good  MTI  operation,  then  p(T)  r 1,  and  consequently 
2 ° 

2(1  -p)a,  <<  1.  This  means  that  in  Eq.  (Al)theterm  exp  \-[p 

2(J.,  2 <l-p2)i 


p‘  = 2 77  (n-m)]  / 
is  negligible  except  when  n = m.  Therefore  Eq.  (Al)  becomes 

on 

.2 


p(p.p') 


277(  1 -pV/2 


exp 


(p-p'Y 

2 3 — 

2o^(\-pZ) 


E exp  - liLi-2-""^^-2"”) 
n=-°°  i (l  + plc/^ 


(A2) 


— 


jBBBE&Utt  PJuQK  MY  *11  Ml) 


jjjgj 
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The  sum  in  Eq.  (A2)  can  be  performed,  approximately,  and  for  p:  1 we  then  find 
from  Eq.  (A2)  that 


P<*.  *'>  - 


* 


exp 


L jsL-r>\\ 

( 4(l-p)ff^J 


where  now  -n  < i/-1  s v and  -is  £ v,  instead  of  their  original  range 
and  -»  £ 4/ ' £ °° . 


i 


i 


(A3) 


s i//s  « 
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Appendix  B 

Clutter  Improvement  Factor  for  Log-Normal  Clutter 


The  ensemble  averaged  residue  power  for  a two-pulse  canceller,  when  clutter- 
only  is  present,  can  be  written  as 


< I r | 2>  = <|Ce^-C'  e^'|2> 


(Bl) 


Upon  recalling  that  C = exp(\),  c'  = exp(x')  we  can  rewrite  Eq.  (Bl)  as 

< | r | 2>  = <e2x+  e2x'  - 2 He  eX+^'+ j . (B2) 


Because  x and  i//  are  gaussian  random  variables  we  can  use  the  well-known  result 
that,  if  0 is  gaussian,  then 

<exp  0)  = expj  <0>  + | < [0  - < 0 > ]2>  } 

If  we  normalize  the  clutter  power  so  that  <exp  2x>  = 1.  and  use  the  aforemen- 

tioned result  we  find  that  Eq.  (B2)  becomes 

<|r|2>  = 2 jl  - He  exp  j\x+x'>  + j <([x-  <X>  + x'  - <x'>  + ^]  [ . 

(B3) 

This  implies  that  <x>  =-  a 2,  where  a 2 = «x  - <x>  )2>  . 

A X 


29 


If  we  now  use  our  previous  assumptions  that  <X*>  = <X>.  (VO  = = 0 • 

2 

X and  are  stationary,  and  the  fact  that  <X>  = -a^'  we  fint*  that  E<4-  (B3) 
becomes 

<|r|2>  = 2 | 1 - exp  [-^x2+  «X'<X>)  <x'  * <x'>  )>  - <^2>  + • 


Upon  recalling  that  < (X  " <X>  ) (X*  ~ < X#>  ^ ~ ^ - O ^ , << ^ ^ 

a^2  p(T),  and  finally  using  our  previous  assumption  in  Eq.  (31)  that  = a , 

we  may  rewrite  Eq.  (114)  as 


<!r|Z>  = 2 jl  - exp  [ - 2 a2  (1  - p)j  j . <B5) 

If  we  remember  that  the  average  gain  in  a two-pulse  canceller  is  2 we  find  that 
the  clutter  improvement  factor  I is 

I = 1 . ( B6) 

1 - exp  (-  2o  ( 1 - p)] 

o 

Finally,  because  p:lwe  have  that  2 a (1  - p)<*  1.  Therefore  we  may  expand 
exp[-2o  tl  - p)]  in  a Taylor  series  to  obtain 

1 r J . (B7) 

2a  (1  -p) 


